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Abstract 

Background: WUSCHEL (WUS)-related liomeobox (WOX) protein family members play important roles in the 
maintenance and proliferation of the stem cell niche in the shoot apical meristem (SAM), root apical meristem 
(RAM), and cambium (CAM). Although the roles of some WOXs in meristematic cell regulation have been well 
studied in annual plants such as Arabidopsis and rice, the expression and function of WOX members in woody plant 
poplars has not been systematically investigated. Here, we present the identification and comprehensive analysis of 
the expression and function of WOXs in Populus tomentosa. 

Results: A genome-wide survey identified 18 WOX encoding sequences in the sequenced genome of Populus 
trichocorpo (PtrWOXs). Phylogenetic and gene structure analysis revealed that these 18 PtrWOXs fall into modern/WUS, 
intermediate, and ancient clades, but that the WOX genes in P. trichocorpo may have expanded differently from the 
WOXqenes in Arobidopsis. In the P. trichocorpo genome, no IVOX members could be closely classified as AtWOXS, 
/\tW0X6, AtW0X7, AtWOXW, and AtWOXM, but there were two copies of 1/1/OX genes that could be classified as PtrWUS, 
PtrWOXI, PtrW0X4, PtrWOXS, PtrWOX8/9, and PtrWOXl 1/12, and three copies of WOX genes that could be classified as 
PtrWOXl and PtrWOXl 3. The use of primers specific for each PtrWOX gene allowed the identification and cloning of 
18 WOX genes from P. tomentoso [PtoWOXs), a poplar species physiologically close to P. trichocorpo. It was found 
that PtoWOXs and PtrWOXs shared very high amino acid sequence identity, and that PtoWOXs could be classified 
identically to PtrWOXs. We revealed that the expression patterns of some PtoWOXs were different to their Arobidopsis 
counterparts. When PtoWOX5o and PtoWOXl 1/1 2o, as well as PtoWUSo ond PtoWOX4o were ectopically expressed in 
transgenic hybrid poplars, the regeneration of adventitious root (AR) was promoted, indicating a functional similarity 
of these four WOXs in AR regeneration. 

Conclusions: This is the first attempt towards a systematical analysis of the function of WOXs in P. tomentoso. A 
diversified expression, yet functional similarity of PtoWOXs in AR regeneration is revealed. Our findings provide useful 
information for further elucidation of the functions and mechanisms of WOXs in the development of poplars. 
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Background 

Homeobox (HB) proteins were first discovered in 
Drosophila, They are a superfamily of transcriptional 
factor proteins containing a conserved 60-amino acid 
homeodomain (HD). HB proteins have been found in 
all eukaryotic organisms tested [1-3]. In plants, a 
great number of HD-containing transcriptional factors 
have been identified in both monocots and dicots [4], 
with KNOTTED 1 being the first identified HD-containing 
protein [5]. The HB protein superfamily is classified into 
six families based on HD sequence, location, association 
with other functional domains, and the protein size and 
structure. These classifications are: homeodomain-leucine 
(HD-Zip), plant homeodomain (PHD)-finger, BELL, 
zinc finger-homeodomain (ZF-HD), WUSCHEL (WUS)- 
related homeobox (WOX), and I<NOTTEDl4il<e-homeobox 
(KNOX) [4]. 

The WOX family of HB proteins is distinguished by the 
phylogenetic relatedness of its homeodomain, and is a plant 
specific HB transcription factor family [6]. The Arabidopsis 
genome encodes at least 15 WOX proteins, and these are 
classified into three clades: [7] a modern/WUS clade 
including WUS and AtWOXl-7; an intermediate clade 
including AtWOXS, 9, 11, and 12; and an ancient clade 
containing AtWOXlO, 13, and 14 [6,8]. Arabidopsis 
WOXs regulate key developmental processes including 
stem cell maintenance in the SAM, RAM, and CAM, 
embryo apical-basal polarity patterning, and lateral organ 
development [8]. AtWUS^ a member of the modern WOX 
clade, is expressed specifically in the organizing center of 
the SAM, and is required for maintaining stem cells via a 
feedback loop with CLV3, a peptide ligand that interacts 
with leucine-rich repeat receptor kinases in neighboring 
cells to restrict the size of the SAM in Arabidopsis [9]. 
AtWOXS y also a member of the modern WOX clade, is 
expressed in root quiescent centre (QC) cells surrounded 
by the stem cells [10]. The AtWOXS protein is essential 
for stem cell maintenance via a negative feedback signal 
provided by CLE40 [11]. AtW0X4 is strongly expressed in 
the CAM, and regulates vascular stem cell maintenance 
[12,13]. AtW0X2 is expressed in zygotes, and is involved 
in regulating cell fate in the apical and basal lineage 
of developing embryos [14]. AtWOX3/PRSl is involved in 
lateral organ development through recruiting organ 
founder cells forming the lateral domain in Arabidopsis 
[15]. AtW0X6 is expressed abundantly in developing 
ovules, and is involved in either ovule patterning or 
differentiation [16]. With regards to stem cell maintenance 
in the SAM and RAM, it is noteworthy that AtWUS and 
AtWOX5 are interchangeable [10]. Furthermore, the 
function of AtWOX3 in lateral organ development can be 
fully complemented by AtWUS [15], and partially comple- 
mented by AtWOX4 [17]. Taken together, it appears there 
is a common mechanism of action among modern WOX 



proteins. It has recently been demonstrated that members 
of the modern WOX clade have evolved, through the 
acquisition of a conserved WUS -box, a repressive activity 
important for leaf blade outgrowth regulation [8]. 

For the intermediate WOXs^ AtWOXS, co-expressed with 
AtW0X2 in the zygote, is required for normal development 
of the pre-embryo [14]. Similar to AtWOXS in Arabidopsis, 
PaWOXS/9 is highly expressed at the early zygotic growth 
stages and the later embryo stages in Picea abies [18]. The 
AtW0X9 gene is involved in maintaining cell division and 
preventing premature differentiation in the Arabidopsis 
SAM [19]. The expression patterns and function of 
AtWOXll and 12 are currently unknown in Arabidopsis, 
but OsWOXll in rice is reported to be required for activat- 
ing shoot-borne root development by directly repressing 
RR2, a cytokinin type- A responsive regulator gene [20]. 

With regards to the most conserved ancient plant WOX 
proteins, the expression and function of AtWOXlO is 
unknown. Although AtWOXlS is expressed ubiquitously it 
is most strongly expressed in developed flowers and young 
siliques, where it is involved in the promotion of replum 
formation during fruit development [21]. AtWOXM is 
thought to prevent premature differentiation of primary 
roots, lateral roots, and floral organs [22]. AtWOXM is pre- 
dominantly expressed in vascular tissues and acts redun- 
dantly with AtW0X4 in vascular cell differentiation [23]. 

Poplars are a widely distributed group of economic 
plants. As a rapid growth species, vegetative growth 
leads to higher production of plant biomass [24], thus 
poplars are considered as feedstocks for bioenergy and 
timber [25,26]. Populus tomentosa is a Chinese native 
poplar species with a high economic value and is widely 
planted in northern China. Because of the key roles that 
WOX proteins play in stem ceU maintenance and lateral 
organ development, WOX proteins are potential targets for 
better and faster growth of P, tomentosa. Here, we identify 
18 WOX encoding genes in P, tomentosa {PtoWOXs), We 
provide a comprehensive analysis of the expression and 
function of the PtoWOXs, We reveal that although 
members of all three WOX clades exist in P, tomentosa, 
PtoWOXs expanded differently from those of Arabidopsis, 
The expression of PtoWOXs was found to differ from their 
Arabidopsis counterparts, with many showing significant 
expression in the roots, and being inducible in the 
regeneration of adventitious roots (ARs). Furthermore, we 
report that the ectopic expression of PtoWOXs from both 
modern and intermediate WOX clades promotes the 
regeneration of ARs. 

Results 

Identification and phylogenetic analysis of WOX genes in 
Populus trichocarpa and R tomentosa 

The P. tomentosa genome has not been sequenced. 
Therefore, to initiate functional analysis of WOXs in 
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P. tomentosa we searched for WOX genes in the 
sequenced genome of P, trichocarpa, a poplar species 
physiologically close to P, tomentosa. We used all 15 
known Arabidopsis WOX protein sequences as queries to 
blast against the P, trichocarpa genome database. Eighteen 
WOX gene sequences were identified in the P, trichocarpa 
genome; this is similar to the number of WOXs found in 
Arabidopsis (15) [7], Oryza sativa (13) [27], Zea mays (18) 
[28], and Vitis vinifera (12) [28]. The putative WOX 
proteins in P, trichocarpa (PtrWOXs) ranged from 171 to 
390 amino acids in length (Table 1). A phylogenetic tree 
of WOX proteins from P, trichocarpa and Arabidopsis 
was generated based on the full WOX amino acid 
sequences (Figure 1, left side). Similar to the HD sequence 
based phylogeny of known plant WOX proteins [29], 
the 18 PtrWOXs fell into three major branches: the 
modern/WUS clade, the intermediate clade, and the 
ancient clade; these were similar to their Arabidopsis 
counterparts (Figure 1, left). The modern/WUS clade 
contained 11 PtrWOXs; these could be divided into five 
subclasses based on their relationship with their Arabidopsis 
counterparts. These were named as PtrWUS, PtrWOXl, 
PtrWOX2, PtrWOX4, and PtrWOX5. The intermediate 
clade consisted of four PtrWOXs; these were classified into 
two subclasses: PtrWOX8/9, and PtrWOXll/12. The 
ancient clade contained only one subclass: PtrWOX13. This 
classification was fiirther confirmed by the exon/intron 

Table 1 Summary of the PtrWOX gene family 



organization of the WOX coding sequences in P. 
trichocarpa and Arabidposis (Figure 1, right). With the 
exception of PtrWOX8/9a and PtrW0Xll/12b, all other 
PtrWOX genes shared a similar intron/exon organization 
in the size and arrangement of introns and exons with 
their corresponding AtWOX counterparts. 

All PtrWOXs contained highly conserved HD residues 
(Additional file 1: Figure SI A). A conserved WUS-box 
domain (TLXLFP) located downstream of the HD domain 
[7] was also found in members of the modern clade, but 
not in the other clades (Additional file 1: Figure SIB). 
In addition, an EAR-like domain was present in the 
C-terminal ends of PtrWUSa, PtrWUSb, PtrWOX5a, 
and PtrWOX5b (Additional file 1: Figure SIC). 

Although the 18 PtrWOXs could be classified into 
modern/WUS clade, intermediate, or ancient clades, it is 
interesting that PtrWOXs may have expanded differently 
in Arabidopsis (Figure 1 left). Two copies of PtrWOXs 
could be classified as WUS (PtrWUSa and b), WOX2 
(PtrWOX2a and b), WOX4 (PtrWOX4a and b), WOX5 
(PtrWOX5a and b), WOX8/9 (PtrWOX8/9a and b), and 
WOXll/12 (PtrWOXl l/12a and b), and three copies of 
PtrWOXs could be classified as WOXl (PtrWOXla, b, 
and c) and WOX13 (PtrWOX13a, b, and c). However, no 
PtrWOXs were closely related to AtWOX3, AtWOX6, 
AtWOX7, AtWOXlO, and AtWOX14 in the phylogenetic 
tree (Figure 1, left). Similarly, close orthologs of AtWOX7, 



Gene name 


Locus 


Genomic position 


Arabidopsis WOXs"" 


Length 


Size (kD) 


PtrWUSa 


Potri.OOSG 114700 


Chr05: 8858336-8859980 


ATIG 17950 


264 


30.004 


PtrWUSb 


Potri.007G012100 


Chr07: 958955-960287 


ATIG 17950 


264 


29.274 


PtrWOXla 


Potri.012G047700 


Chrl2: 4446344-4449088 


AT3G18010 


387 


43.851 


PtrWOXlb 


Potri.015G039100 


Chrl 5: 3566414-3569217 


AT3G18010 


374 


42.326 


PtrWOXlc 


Potri.OlOGl 11400 


ChrlO: 13009243-13010950 


AT3G18010 


316 


35.699 


PtrWOXla 


Potri.001G237900 


ChrO 1:24918823-24920230 


AT5G59340 


245 


27.877 


PtrWOXlb 


Potri.009G029200 


Chr09: 3976470-3977823 


AT5G59340 


246 


27.601 


PtrW0X4a 


Potri.002G124100 


Chr02: 9309765-9311320 


AT1G46480 


213 


24.417 


PtrW0X4b 


Potri.014G025300 


Chrl4: 2169818-2171379 


AT1G46480 


213 


24.463 


PtrWOXSa 


Potri.008G065400 


Chr08: 3963469-3964104 


AT3G 11160 


181 


21.009 


PtrWOXSb 


Potri.OlOGl 921 00 


ChrlO: 18688337-18689050 


AT3G11160 


171 


19.615 


PtrWOX8/9a 


Potri.004G051600 


Chr04: 4008260-4010297 


AT1G33880 


390 


43.072 


PtrWOXmb 


Potri.011G061400 


Chrl 1: 5466036-5468417 


AT1G33880 


377 


41.601 


PtrWOXl 1/1 la 


Potri.013G066900 


Chrl 3: 5239050-5240692 


AT3G03660 


255 


28.086 








AT5G17810 






PtrWOXl 1/1 lb 


Potri.019G040800 


Chrl 9: 4693429-4695238 


AT3G03660 


226 


25.001 








At5G17810 






PtrWOXlSa 


Potri.005G101800 


Chr05: 7788135-7790958 


AT4G35550 


248 


27.984 


PtrWOXlSb 


Potri.005G252800 


Chr05: 25454469-25456507 


AT4G35550 


216 


24.639 


PtrWOXlSc 


Potri.002G008800 


Chr02: 484202-486345 


AT4G35550 


215 


24.479 



^Arabidopsis l/l/OX genes that are closely classified with PtrWOXs according to the phylogenetic analysis of protein sequences. 
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PtrWUSa 
PtrWUSb 
AtWUS 
PtrWOXSa 
PtrWOXSb 
AtWOXS 
MSNOXl 
PtrWOXIa 
PtrWOXIb 
PtrWOXIc 
AtWOXI 

Atwoxe 

PtrW0X4a 
PtrW0X4b 
AtW0X4 
AtWOXS 
PtrW0X2a 
PtrW0X2b 
AtW0X2 
AtWOXS 
AtW0X9 
PtrWOX8/9a ^ 
PtrWOX8/9b = 
PtrW0X11/12a^ 
PtrW0X11/12bL- 
AtW0X11 
AtWOXI 2 
PtrW0X13a = 
AtWOXI 3 
PtrWOXISb ^ 
PtrW0X13c 
AtWOXI 4 
AtWOXI 0 



500bp 



Figure 1 Phylogenetic tree of WOX proteins and gene structure 
of corresponding WOX genes from Populus trlchocarpa and 
Arabidopsis thaliana. Deduced full-length amino acid sequences 
were aligned using ClustalX 2.0. The phylogenetic tree was constructed 
using phyML and the maximum likelihood method. The three clades 
(ancient, intermediate, and modern/WUS) are indicated as different 
grayscale boxes in the left panel. Empty boxes and black lines in the 
right panel represent exons and introns respectively. Support values 
are shown on selected branches. 



AtWOXlO, AtWOX12, and AtWOX14 were missing from 
Vitis vinifera [28], Picea abies [30], and other woody plant 
species [28]. Details of these PtrWOX encoding genes, 
including their corresponding Arabidopsis counterparts, 
and amino acid lengths are listed in Table 1. 

Using primers specific to each PtrWOX gene (Additional 
file 2: Table SI), we identified and cloned 18 correspond- 
ing WOX genes from P, tomentosa {PtoWOXs), We found 
that PtoWOXs and PtrWOXs shared very high amino 
acid sequence similarities (86-99%) and CDS sequence 
identities (Additional file 2: Table S2). Phylogenetic 
analysis revealed that the PtoWOXs could be classified 
identically to PtrWOXs (Figure 2). 

Chromosomal location and gene duplication of WOX 
genes in R trichocarpa 

To better understand how WOX genes expanded in 
poplars, we took advantage of the sequenced genome 
of P. trichocarpa and mapped the chromosomal location 
of the 18 PtrWOX loci. These 18 PtrWOX sequences were 



56 



53 



99 



100 



98| PtoWUSa 



100 



100 



PtrWUSa 
PtoWUSb 
75L PtrWUSb 
96 p PtoWOX5a 
>- PtrW0X5a 

I- PtoWOX5b 

99I- PtrW0X5b 
PtoWOXIc 



L PtrWOXIc 
PtoWOXIa 
PtrWOXIa 
PtoWOXIb 



100 

96 1 PtrWOXIb 



77rPtoWOX4a 
P- PtrW0X4a 
lOoL PtoWOX4b 
89l PtrW0X4b 



100 



lOOr PtoWOX2a 



' PtrW0X2a 
PtoWOX2b 
99' PtrW0X2b 

100|PtoWOX13b 



99 



100 



PtrWOXISb 
iPtoWOXISc 



100lptrWOX13c 
iPtoWOXISa 



100"PtrWOX13a 
99 | PtoWOX8/9a 

IO0J I- PtrWOX8/9a 

[1 PtoWOX8/9b 
60' PtrWOX8/9b 
99 59rPtoWOX11/12a 

[L PtrWOX1 1 /1 2a 

T00|_r PtoWOX11/12b 
751 PtrW0X11/12b 

Figure 2 Phylogenetic tree of WOX family proteins from P. 

trichocarpa and P. tomentosa. Multiple alignment of full-length 
WOX protein sequences from the two poplar species was executed 
using Clustal X2.0, and a phylogenetic tree constructed using MEGA 
4.0 and the neighbor-joining (NJ) method with 1000 bootstrap replicates. 
Bootstrap support is indicated at each node. 



0.1 



unevenly distributed among the 14 chromosomes, with 
the exception of chromosomes III, VI, XIII, XVI, and 
XVIII (Additional file 1: Figure S2). The chromosomal 
duplication map of P, trichocarpa was generated pre- 
viously [31]. We identified five pairs of PtrWOX 
genes, PtrWOX la-PtrWOX lb, PtrWOX4a-PtrWOX4b, 
PtrWOXSa-PtrWOXSh PtrWOXll/12a-PtrWOXll/12b, and 
PtrWOXlSb-PtrWOXlSc in the duplicated segments of 
the P, trichocarpa genome. This suggests that expansion 
of the PtrWOX gene family was, at least partially, caused 
by chromosomal duplication events during the evolution 
of poplar species. 

Nuclear localization of PtoWOXs 

As a first step towards understanding the functions 
of WOXs in P, tomentosa y we examined the subcellu- 
lar localization of PtoWUSa, PtoWOX4a, PtoWOXSa, 
PtoWOXll/12a, and PtoWOXlSc members of selected 
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PtoWOX subclasses. As transcription factors, WOXs 
should be targeted to the nucleus. To test this, we fused 
YFP at N-termini of PtoWUSa, PtoWOX4a, PtoWOXSa, 
and PtoWOXll/12a and the C-terminus of PtoWOXlSc. 
In a transient expression assay of Nicotiana benthamiana 
leaf epidermal cells all fusion proteins were found in the 
nucleus (Figure 3). 

Diversified expression patterns of PtoWOX genes 

To gain insights into possible developmental and 
physiological functions of WOX proteins in P, tomentosa, 
the expression pattern of WOX genes in P, tomentosa Carr. 
were analyzed by semi-quantitative RT-PCR (Figure 4, 
Table 2). Two PtoWOX2s and two PtoWOXll/12 s were 
mainly expressed in the roots. In contrast, PtoWUSb, and 
three PtoWOXls were nearly absent from roots. Further to 
this, PtoWOXla and lb were only expressed in leaves and 
PtoWOXlc only in the stem. PtoWOXSa and PtoWOX8/9a 



A 












35S::YFP 
-PtoWUSa 


briqhtfield 




B 


♦ 




• 






35S::YFP 
-PtoWOX4a 


brightfield 


✓ * 

Merge — 


C 






35S::YFP 
-PtoWOXSa 




Merge 


D 






• 


■ 


e 


35S::YFP- 
PtoWOX11/i: 






E 






m 




E 


35S::PtoWO> 
-YFP 


^^%rightfield 





Figure 3 Subcellular localization of five PtoWOX proteins. The 

nuclear localization of PtoWUSa (A), PtoWOX4a (B), PtoWOXSa (C), 
PtoWOXl 1/1 2a (D), and PtoWOXlSc (E) in a tobacco lower 
epidermal cell. The YFP channel, the bright field, and the merged 
images are shown on the left, middle, and right panels respectively. 



were highly expressed in both roots and leaves, while 
PtoWOX8/9b was expressed in both roots and stem. 
It is noteworthy that PtoWUSa, PtoWUSb, PtoWOX4a, 
and PtoWOX4b were strongly expressed in the CAM zone 
of poplar stems. In Arabidopsis, AtWUS is only expressed 
in the SAM [32]. In addition to the enrichment of 
PtoWUSa and PtoWUSb in the CAM zone, PtoWUSa 
was also expressed, albeit weakly, in roots, leaves, and 
stems. Similarly, AtWOXS is restricted to the RAM [10], 
while PtoWOXSb was ubiquitously expressed. The only 
subclass of the ancient WOX clade in poplars, PtoWOXlSs 
{PtoWOXlSa, b, c), also exhibited a ubiquitous expression 
pattern. In a separate research project conducted in our 
laboratory, we generated a set of RNA-seq data for 
the hybrid poplar Populus alba X Populus glandulosa. 
Quantification of expression of the corresponding 18 
WOX genes in P, alba X P, glandulosa {PagWOXs; 

c 



PtoWUSa 

PtoWUSb ^^^^^^^^^1 

PtoWOXla ^BSHHI^^I 
PtoWOXIb ^^^^^^^^^B 

PtoWOX2a BHHHm^H 
PtoWOX2b ^^^^^B^^^l 
PtoWOXla ^l^^^l^^^l 

PtoWOX4a ^^^^^^^^^3 
PtoWOX4b ^^^^^^^^^g 

PtoWOX5a ^^^^l^^^^l 

PtoWOXSb ^^^^^^^^^1 

PtoWOX8/9a ^^^9^H^^I 
PtoWOX8/9b ^^^^^^^^^H 

PtoWOX 11 /12a ^^i^^mm 

PtoWOX11/12b BI^BI^^^^I 
PtoWOX13a ^B^^^^^^H 
PtoWOXISb ^^^^^^^^^B 
PtoWOXlSc ^^^^^^^^^3 
PtoActin P^^^^^^^^^^3 

Figure 4 Semi-quantitative analysis of Pto I/I/OX genes. The 

expression of PtoWOXs in roots, leaf, stem, CAM zone, and immature 
xylem from P. tomentoso Can. The amplification cycle used was 32 
for PtoWOX4o, 4b, and 13c, 35 for PtoWUSo, PtoWOXlo, 13o and 13b, 
and / 1/120 and / l/12b; 40 for all other WOXs, and 30 for the 
PtoActin reference control. 

V J 
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Table 2 Summary of the expression pattern of AtWOXs and PtoWOXs 


PtoWOXs 


Expression pattern 


AtPWOXs 


Expression pattern 


Ref 


PtoWUSo 


Root, leaf, SAM, cambium, xylem 


AtWUS 


SAM 


[32] 


PtoWUSb 


Cambium, xylem 








PtoWOXlo 
PtoWOXlb 


Leaf 
Leaf 


AtWOXl 


Between the adaxial and abaxial 
domains of leaf 


[33] 


PtoWOXlc 


Stem 








PtoWOX2o 


Root 


AtVvUX2 


Apical embryo domain 


[34] 


PtoWOX2b 


Root 








No closely classified counterparts 


AtVvUAj 


flower primordia, floral organ primordia, 
and young leaf primordia. 




PtoWOX4a 
PtoWOX4b 


Root, Leaf, stem, cambium, xylem 


AtW0X4 


Cambium, trichomes, stomata, phloem, 
pericycle 


[12] 


PtoWOXSo 


RAM, leaf 


AtWOXS 


RAM 


[10] 


PtoWOXSb 


Root, leaf, stem 








No closely classified counterparts 


AtW0X6 


Differentiating primordia and developing ovules 


[16] 


No closely classified counterparts 


AtWOX/ 


unknown 




PtoWOX8/9a 


Root, leaf 


AtWOXS 


zygote, proembryo and embryo 


[14] 


PtoWOX8/9b 


Root, stem 


AtW0X9 


Developing embryos, proliferating tissues, SAM, 
leaf primordia, floral meristems, epidermal layer 
of the placenta and growing septum 


[19] 


No closely classified counterparts 


AtWOXlO 


unknow 




PtoWOXl 1/120 


Root 


OsWOXl 1 


Cell division regions in roots and shoots. 


[20] 


PtoWOXl 1/1 2b 




AtWOXl 1, 12 


unknown 




PtoWOXlSo 
PtoWOXlSb 


Root, leaf, SAM, cambium, xylem 


AtWOXlS 


Root tip, emerging lateral roots, root, SAM 
and leave vasculature, gynoecia 


[21] 


PtoWOXlSc 










No closely classified counterparts 


AtW0X14 


Vascular tissue 


[23] 



The l/l/OX expression patterns in other plant species were published In the corresponding references listed. Tissue names In Italic Indicate strong expression. 



Additional file 1: Figure S3, Additional file 2: Table S3) 
indicated that the expression pattern of 18 PagWOX genes 
in the hybrid poplar was similar to that of the 18 PtoWOX 
genes in P. tomentosa, though there was a discrepancy 
between PtoWUSb and PagWUSb, 

Promoter::GUS based analysis of expression of PtoWOX 
genes 

RT-PCR is a rapid way to examine tissue specificity of 
gene expression, but it does not provide a good resolution. 
Since we are interested in the process of the root growth 
and wood formation in Populus tomentosa, a root- 
specific PtoWOXll/12a, a root-enriched PtoWOXSa, and 
two generally expressed but CAM enriched PtoWOXs 
(PtoWUSa and PtoWOX4a) were selected for promoter:: 
GUS assays to examine their expressions in detail Due to 
ease of genetic manipulation, the hybrid poplar P, alba 
X P, glandulosa, a close relative of P. tomentosa, was 
chosen for gene transformation. In Arabidopsis, AtWOXS 
is expressed specifically in the QC cells in the root tips 



[10]. Semi-quantitative RT-PCR results revealed that in 
addition to being strongly expressed in the roots 
PtoWOXSa was weakly expressed in the leaves (Figure 4). 
The promoter ::GUS assays indicated that PtoWOXSa was 
mainly expressed in a small area behind the cap region of 
ARs in the Pptowoxsa''^^^ transformed hybrid poplar 
(Figure 5B, see arrow). The PtoWOXSa promoter had 
a weak activity in the leaves (Figures 4 and 5A). 
Semi-quantitative RT-PCR identified PtoWOXll/12a as a 
root specific WOX gene in P, tomentosa (Figure 4). A 
GUS assay confirmed it was only expressed in the roots, 
especially in a small region behind the AR cap (Figure 5C 
and D, see arrow). Based on The semi-quantitative 
expression analysis revealed PtoWUSa and PtoPWOX4a 
were strongly expressed in the CAM zone, but were 
detectable in other examined tissues (Figure 4). GUS 
staining (Figure 5E-J) revealed that PtoWUSa and 
PtoPWOX4a were both strongly expressed in the 
CAM zone of stems (Figure 5G and J); they were also 
weakly expressed in roots (Figure 5F and I) as GUS 
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PtoWOX4a- 



Figure 5 GUS-staining assays of four PtI/l/OX promoters. (A-B) GUS expression in (A) 2-weel<-old sapling and (B) root of Ppfowoxja-G^j'S-expressing 
poplar plants. The insert in A shows GUS expression in the leaf, and that the GUS expression in B is mainly focused on a small area immediately behind 
the root cap in AR tips. (C-D) GUS expression in (C) 2-week-old sapling and (D) root tip of Ppfo^ox////2o"GL'S-expressing poplar plants. GUS is expressed 
in a small area immediately behind the root cap in AR tips, and the insert in C shows no GUS expression in leaf (E-G) GUS expression in (E) 2-week-old 
sapling, (F) root including a small area behind the root cap indicated by an arrow, and (G) the CAM zone of Pptowuso-G^'S-expressing poplar plants. The 
insert in E shows GUS expression in the SAM area in a halved shoot apical part. (H-J) GUS expression in (H) 2-week-old sapling, (I) root tip and (J) stem 
of ^Ptoi/^ox4a"G6'S-expressing poplar plants. The insert in H shows no GUS expression in the SAM area in a halved shoot apical part. The bar represents 
1 cm in A, C, E, and H; 1 mm in D; 0.5 mm in B, F, and I; and 0.1 mm in G, J, and the inset of E. 



staining of roots of transgenic lines expressing either Pp^o- 
wusa-G^^ or Pptowox4a''G^S occurred after a prolonged 
period. Interestingly, in the Pptowusa--^^^ lines, a small 
region behind the root cap could be stained (Figure 5F, 
arrow), though activity of Pptowusa-'^^^^ not limited to 
this small region in roots. In Pptowox4a''G^^ transgenic 
lines, the small region behind the root cap region was not 
stained, but the root cap region of ARs could be stained 
(Figure 51, arrowhead). 



Dynamic expression of PtoWOXs during the regeneration 
of adventitious shoots and roots in R tomentosa 

The regeneration of adventitious shoot (AS) and AR using 
leaf or stem explants is one of the best tools for rapid 
propagation and genetic manipulation of poplars. These 
regeneration processes involve the re- establishment of 
SAM and RAM, which leads to the final differenti- 
ation of leaves and roots. Therefore regeneration of 
AS and AR is an excellent system for studying stem 
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cell initiation, cell fate determination, and hormonal 
signaling [36]. AS regeneration is composed of several 
successive stages, including the Bl pre-induction stage 
(Additional file 1: Figure S4A), the B2 callus induction 
stage (Additional file 1: Figure S4B), the B3 callus expan- 
sion stage (Additional file 1: Figure S4C), the B4 callus 
transition stage (Additional file 1: Figure S4D), the B5 AS 
formation stage (Additional file 1: Figure S4E), and 
the B6 AS growth stage (Additional file 1: Figure S4F). 
To investigate whether PtoWOXs play roles in the 
regeneration of AS we examined the dynamics of PtoWOX 
expression using quantitative real-time PGR. Among 
the 18 PtoWOXs, 10 had significant changes in expression 
during AS regeneration (Figure 6A). The expression 
levels of PtoWUSa, PtoWOX4a, PtoWOX4b, PtoWOXSa, 
and PtoWOXSb changed over twenty-fold during the 



regeneration of AS. They had different expression dynam- 
ics with maximum expression of PtoWOX4a/b seen at 
earlier stages (B2/B3), and PtoWOXSa/b at later stages 
(B4) (Figure 6). All three PtoWOXlSs were expressed at a 
relatively low level in AS regeneration but were slightly 
induced during the AS regeneration process (Figure 6A). 
In contrast, expression of PtoWOXla and lb was slightly 
down-regulated in AS regeneration. This down-regulation 
of PtoWOXls may be a technical artifact as both 
PtoWOXls were expressed specifically in leaves, thus 
their expression in the Bl stage may be still at a high level. 
Nevertheless, we view this as an indication that PtoWOXls 
are not involved in SAM and AS regeneration. Finally, 
transcription of PtoWOX2a, PtoWOX2b, PtoWOXll/12a, 
and PtoWOXll/12b, whose expression was mainly in the 
roots, was unchanged during AS regeneration. 
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Figure 6 Quantitative PCR analysis of dynamic expression of PtoWOXs in the regeneration of AS and AR. (A) Dynamic changes of PtoWOXs in AS 
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AR regeneration is composed of a Rl pre-induction stage 
(Additional file 1: Figure S4G), R2 callus formation stage 
(Additional file 1: Figure S4H), R3 AR emergence stage 
(Additional file 1: Figure S4I), and a R4 AR elongation 
stage (Additional file 1: Figure S4J). We examined the 
dynamics of PtoWOXs expression to gain insights into 
their possible roles in these processes. Expression of 14 out 
of 18 PtoWOXs was detected and changed during the 
AR regeneration process (Figure 6B). Among these, 
expression of four WOX genes: PtoWOXSa, PtoWOXSb, 
PtoWOXll/12a, and PtoWOXll/12b changed significantly 
during AR regeneration (Figure 6B). The expression of 
PtoWOXll/12a was over 80-fold stronger and maintained 
a high level in the R2 though R4 stages. PtoWOXll/12b 
was induced at the R3 and R4 stages, but not at the Rl or 
R2 stage. It is interesting that although expression of 
PtoWOXSa and PtoWOXSb was induced at the R2 
stage, their expression gradually decreased in the following 
stages. Similarly, root specific WOXs were not induced in 
AS regeneration; the expression levels of PtoWOXla and 
PtoWOXlb, two leaf-specific WOX genes in P. tomentosa 
(Figure 4), were not changed. 

Ectopic expression of PtoWOXSa and PtoWOX11/12a as 
along with PtoWUSa and PtoWOX4a promotes AR 
regeneration 

Because of the high expression of PtoWOXSs, PtoWOXll/ 
12 5, PtoWUSa, and PtoWOX4a in the roots, and the strong 
induction of PtoWOXSs and PtoWOXll/ 12 s in the 
regeneration of AR, we were interested in the roles 
they play in the AR regeneration process. We used a 
35S promoter to ectopically express PtoWOXSa and 
PtoWOXll/12a along with PtoWUSa and PtoWOX4a in 
the hybrid poplar {P. alba X P. glandulosa), AR regenerated 
in stem segments of three independent transgenic 
plants were quantified at the young sapling stage. All 
four ectopically expressed WOX genes promoted AR 
regeneration in the transgenic hybrid poplar lines (Figure 7). 
It appeared that PtoWOXSa and PtoWOXll/12a had 
a better ability than PtoWUSa and PtoWOX4a to promote 
the regeneration of AR. This is consistent with the strong 
induction of PtoWOXSa and PtoWOXll/ 12a in the 
regeneration of AR. Although the expression of PtoWUSa 
and PtoWOX4a was not obviously induced (Figure 6B), 
their effect on the regeneration of ARs when ectopically 
expressed indicated that the four PtoWOX genes may have 
overlapping functions. 

Discussion 

Diversification of WOX genes In P. trichocarpa and P. 
tomentosa 

We identified 18 WOX encoding genes in both P. 
trichocarpa and P, tomentosa. Although these poplar 
WOX genes could be categorized into modern/WUS, 



intermediate, and ancient clades, it is interesting that 
no WOX genes from P, trichocarpa and P, tomentosa are 
classified together with AtWOXS, AtW0X6, AtWOXV, 
AtWOXlO, or AtW0X14, Furthermore, many of the 
WOX genes in P, trichocarpa had sister copies owing 
to chromosomal duplication events. 

While the roles oi AtW0X6 and AtWOX? are not well 
defined, AtWOXS is expressed in leaves and involved in 
leaf blade outgrowth [15,37]. Moreover, the function of 
the AtWOX3 protein can be fully complemented by 
AtWUS [15], and partially complemented by AtWOX4 
[17]. It is possible that the WOX3-mediated function of 
leaf blade outgrowth regulation is compensated by other 
modern WOX members in poplars. In this regard, we 
note that the PtoWOXla gene and PtoWOXlb gene are 
specifically expressed in P, tomentosa leaves (Figure 4). 

Arabidopsis contains three ancient WOX genes, AtWOXlOy 
AtWOXlS, and AtWOXM. P. trichocarpa and P, tomentosa 
also have three ancient WOX genes, PtrWOXlSa, b, and c; 
these are very similar in sequence, with PtrWOXlSb and 
PtrWOXlSc being sister pairs. Grape also has three ancient 
WOX genes, VvWOXlSA, B, and C. [28]. There appears to 
be functional diversification in AtWOXlO, AtWOXlS, and 
AtWOXM; the AtWOXlS gene is involved in replum for- 
mation during fruit development [21], AtWOXM acts re- 
dundantly with AtW0X4 in vascular cell division [23], 
while the role of AtWOX 10 is unknown. The three ancient 
WOX genes of P. trichocarpa and P. tomentosa are ubi- 
quitously expressed in all examined tissue in poplars 
(Figure 4). It is possible that the diversified functions 
of ancient AtWOXlO, 13, and i4 in Arabidopsis may 
be maintained by three ancient WOX genes, namely 
WOXlSa, by and c in poplars. 

Expression of PtoWOX genes 

Protein function can be specified in terms of temporally 
and spatially regulated gene expression. We note that 
some PtoWOXs genes are expressed differently from 
their Arabidopsis counterparts (Table 2). AtWUS is 
expressed specifically in the rib-meristem (RM) cells 
beneath the central zone of the SAM, and it maintains the 
stem cell population via a regulatory loop with CLAVATA 
genes in Arabidopsis [9,32]. In poplar, the expression of 
PtoWUSa is not limited to the SAM and is strongly 
induced in the regeneration of AS; it is also detected 
in roots, stem and leaves (Figure 4, Figure 5E-G), and 
particularly in the CAM zone of stems. A major difference 
between Arabidopsis and poplar is that the latter is a 
perennial plant with annual differentiation of vascular 
tissues from CAM. It seems that the PtoWUSa gene 
may also be involved in vascular tissue differentiation. 
In Arabidopsis y the AtW0X4 gene is required for vascular 
differentiation from CAM cells [12,38]. Both PtoWOX4a 
and PtoWOX4b genes are strongly expressed in CAM cells 
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Figure 7 Ectopic over-expression of four PtoWOXs promotes regeneration of AR. (A-B, E, and F) Regenerated ARs in transgenic poplar 
plants ectopically expressing (A) PtoWUSa, (B) PtoWOX4a, (E) PtoWOXSa, and (F) PtoWOXl 1/1 2a. Wild type and over-expressing lines are on the 
left and right of each panel, respectively. (C-D, G, and H) Statistical result of ARs in different transgenic lines. No. of AR: adventitious root numbers; bars 
in A, B, E, and F represent 1 cm. '''' in C, D, G, and H indicates significant difference using LSD test at P <0.01 . 



and are strongly induced in AS regeneration, but not in 
AR regeneration (Figure 6). Therefore, differently from 
PtoWUSa, the function of the PtoWOX4 gene in CAM 
cells may be maintained in poplars. It will be interesting 
to test whether PtoWUSa and PtoWOX4 act differently in 
vascular tissue differentiation in poplars. 

AtWOXS is a QC cell specific gene in Arabidopsis [11]. It 
seems that functional expansion also occurs for PtoWOXSs. 
In poplars, the PtoWOXSa gene is strongly expressed in 
a small region that resembles QC cells, but is also detect- 
able in young leaves (Figure 4). Furthermore, PtoWOXSa 



expression is strongly induced in the regeneration of AR 
as well as in the AS. Therefore, we speculate that the 
function of PtoWOXSa and PtoWOXSb may have 
expanded from roots to leaves in poplars. 

Although the expression pattern of AtWOXll in 
Arabidopsis has not been reported, expression of OsWOXll 
in rice is detected in cell division regions of both roots 
and shoots [20]. We note that expression of both 
PtoWOXl 1/12 s is restricted to a small area behind 
the root cap region that resembles the root QC in AR tips, 
and probably the pericycle (Figure 5D). Thus, the action 
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of PtoWOXll/12 s may be restricted to root development 
in poplars. During the regeneration of AR, the expression 
levels of PtoWOXll/12 s and PtoWOXSs are all strongly 
induced; however, it is interesting that the expression of 
PtoWOXll/12a is induced at the R2 stage, while the 
expression of PtoWOXll/12b is induced at the R3 stage. 
It is likely that PtoWOXll/12a acts differently from 
PtoWOXll/12b in the regeneration of AR. 

Four PtoWOX proteins can perform similar functions in 
promoting AR regeneration 

A key function of WOX proteins is to maintain the stem 
cell population in different tissues and organs [9,11,12,39]. 
Recently, it has been demonstrated [8] by using a laml 
complementation screen of WOX proteins in N. sylvestris 
that modern members of WOX proteins may have 
acquired additional functions during evolution that lead to 
functional specificity of this clade in the regulation of 
leaf blade outgrowth. In this study, we found that when 
ectopically expressed, three members of the modern 
WOXs (PtoWOX4a, PtoWOXSa, and PtoWUSa) and a 
member of the intermediate WOXs (PtoWOXll/12a) 
promote AR regeneration (Figure 7). In our promoter:: 
GUS based expression analysis, the PtoWUSa, PtoWOXSa, 
and PtoWOXll/12a genes are expressed in a small region 
in the root tips resembling QC cells, the PtoWOX4a gene 
however, was not expressed in restricted areas in the root 
tips. Therefore, it is likely that members of the modern 
and intermediate WOXs can perform similar functions in 
the regeneration of AR in poplars, even if the gene is not 
normally expressed in root tip cells. 

There is a clear divergence in the amino acid sequences 
of the 15 Arabidopsis and 18 poplar WOX proteins, with 
only modern WOX proteins possessing a WUS-box 
domain. Lin et al [8] demonstrated that the WUS-box is 
required for the regulation of leaf blade outgrowth. Many 
different WOXs act in different cells and tissues to 
maintain the stem cell population [9,11,12,39], thus a 
common action mechanism for WOX proteins in 
maintaining the stem cell niche has been proposed 
[8]. Our results add support to this viewpoint. The 
regeneration of AR relies on a balance between cell 
differentiation and renewal of stem cells in the RAM 
[11]. It is possible that PtoWOXSa and PtoWOXll/12a, 
whose genes are expressed and strongly induced in the 
regeneration of AR, are major players in maintaining 
the stem cell niche in root tips. However, the function of 
PtoWOXSa and PtoWOXll/12a in maintaining the 
stem cell niche in root tips may be partially overlapped by 
PtoWUSa or PtoWOX4a. When the latter two PtoWOX 
genes are ectopically expressed they promote the regener- 
ation of AR. This functional overlapping may be attributed 
to a common mechanism for maintaining the stem cell 
niche. It is possible that a WOX domain other than the 



WUS-box is responsible for this common action. It 
will be interesting to examine which WOX domains 
have evolutionary significance in the function of WOXs in 
stem cell maintenance. 

Conclusions 

This study represents a step forward in our understanding 
of the functions and mechanisms of 18 WOXs found in 
poplar, in particular, their possible roles in vascular as well 
as root development of poplars. Based on the expression 
of PtoWOX genes revealed in this study, it is attempted to 
suggest that, in poplars, PtoWOX4s and PtoWUSs are 
involved in vascular development, while PtoWOXSs and 
ProWOXll/12s are major players in root development. 
Despite these possible functional specificities, it appears 
that there is a common action mechanism for different 
PtoWOXs in maintaining different stem cell niches. 

Methods 

Bioinformatic analysis 

The P, trichocarpa genome (release 3.0, http://www. 
phytozome.net/poplar) was blasted using AtWOXs pro- 
tein sequences as queries with NCBI BLASTP. Obtained 
sequences were used as secondary queries to re-blast the P, 
trichocarpa genome. After removing redundant sequences, 
multiple alignments of full length sequences of AtWOX 
and PtrWOX proteins were performed using the Clustal 
X2.0 program [40]. Maximum likelihood (ML) phylogenetic 
trees were constructed using PhyML (v3.0) with JTT 
amino acid substitution model, 1000 bootstrap replicates, 
estimated proportions of invariable sites, estimated gamma 
distribution parameters, and an optimized starting BIONJ 
tree [41,42]. A multiple alignment of full-length WOX 
protein sequences from P, trichocarpa and P, tomentosa 
was executed using Clustal X2.0 [40], and a phylogenetic 
tree constructed using MEGA 4.0 by the neighbor-joining 
(NJ) method with 1000 bootstrap replicates. The trees 
have been submitted to Treebase under study number 
1S612 (http://treebase.org/treebase-web/search/study/ 
summary.html?id=lS612). Accession numbers used in 
this study are listed in Additional file 2: Table S4. Exon 
and intron structures of individual PtrWOXs were illus- 
trated using the Gene Structure Display Server (GSDS, 
http://gsds.cbi.pku.edu.cn/) [43] by aligning the cDNA se- 
quences with the corresponding genomic DNA sequences 
from www.phytozome.net. All 18 PtrWOX genes were 
mapped to P, trichocarpa chromosomes. Whole-genome 
duplication analyses were accomplished as described in 
Tuskan et al [30]. The MEME program (version 4.3.0, 
http://meme.sdsc.edu) [44] was used for elucidation of 
motifs in HDs. MEME was run locally with the following 
parameters: number of repetitions - any; maximum num- 
ber of motifs - 20; and the optimum motif widths were 
constrained from 6 to 21 residues. 
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Plasmids and constructs 

The coding sequences of all PtoWOXs except for PtoWUSb 
(submitted to NCBI by another group) were amplified 
from the cDNA of P. tomentosa, and cloned into 
pDNOR222.1 (Life technologies, Carlsbad, California, 
U.S.) (to produce pENTRs) for sequencing. PtoWUSa, 
PtoWOX4a, PtoWOXSa, and PtoWOXll/12a were sub- 
cloned into pMDC32 to produce 35S::PtoWUSa, 35S:: 
PtoWOX4a, 3SS::PtoWOXSa, and 3SS::PtoWOXll/12a 
constructs. The resultant pENTR constructs were then 
recombined into pEarleyGatel04 (ABRC stock DB3-686) 
to produce 3SS::YFP-PtoWUSa, 3SS::YFP-PtoWOX4a, 
3SS::YFP-PtoWOXSa, and 3SS::YFP-PtoWOXll/12a con- 
structs using the Gateway cloning system (Life tech- 
nologies, Carlsbad, California, U.S.). PtoWOX13c without 
the stop codon was amplified and subcloned into pEarley- 
GatelOl (ABRC stock DB3-683) to produce 35S::Pto 
W0X13C-YFP, The primer sequences used for amplifica- 
tion of PtoWOXs are listed in Additional file 2: Table SI. 
5'UTR fragments, 2-^3 kb in size, of PtoWUSa, Pto- 
W0X4a, PtoWOXSa, and PtoWOXll/12a were amplified 
from the genomic DNA of P. tomentosa Carr. Primer se- 
quences and promoter lengths are listed in Additional file 2: 
Table SI. The amplified promoter fragments were cloned 
into pDNOR222.1, and then subcloned into pMDC164 to 

produce PptoWUSa-'GUS, PptoWOX4a''GUS, PptoWOXSa-'GUS, 

and Pptowoxii/i2a''GUS constructs using the gateway clon- 
ing system (Life technologies, Carlsbad, California, U.S.). At 
least three independent lines were used for analysis. 

Plant cultivation and transformation 

Tobacco plants {Nicotiana benthamiana) used for transient 
expression were grown on soil under an 8/16 h (day/night) 
photoperiod at 20°C. All constructs were introduced 
into Agrobacterium tumefaciens strain GV3101 by 
electroporation. A single A, tumefaciens colony con- 
taining 3SS::YFP-PtoWUSa, 3SS::YFP-PtoWOX4a, 3SS:: 
YFP-PtoWOXSa, 3SS::YFP-PtoWOXll/12a, or 3SS:: 
PtoWOX13c-YFP was used to inoculate 2 mL of YEP 
medium (per liter: 10 g tryptone, 10 g yeast extract, 
5 g NaCl, pH 7.0), supplemented with 50 mg/L kanamycin, 
10 mg/L gentamycin, and 34 mg/L rifampicin. Bacterial 
cultures were incubated at 28°C with agitation until 
OD6oo = 0-5- 0.5 mL of culture was transferred into an 
Eppendorf tube, and the bacteria pelleted by centrifugation 
at 2000 X ^ for 5 min in a microcentrifuge at room 
temperature. The pellet was washed twice with 0.5 mL of 
infiltration buffer (10 mM MgCl2, 150 (iM acetosyrin- 
gone) and resuspended in 0.5 mL of the same buffer. The 
inoculum concentration of 3SS::YFP-PtoWUSa, 35S::YFP- 
PtoWOX4a, 3SS::YF-PtoWOX5a, 3SS::YFP-PtoWOXll/12a, 
and 3SS::PtoWOX13c-YFP (ODeoo = 0.1) was adjusted by 
diluting the bacterial suspension with the infiltration buf- 
fer. The inoculum was delivered to tobacco lower 



epidermal leaf cells by gentle pressure infiltration using a 
1-mL syringe without a needle. Following infiltration, plants 
were incubated under dark conditions for 6 h, and then 
grown under normal conditions, as described above. To- 
bacco leaves were analyzed 3 days after infiltration. 

Hybrid poplar {P. alba X P. glandulosa) clone 84K 
used for transformation were kept at 23-25°C under a 
16/8 h (day/night) photoperiod, with light intensity of 
50 (iM m"^ s"^ provided by cool white fluorescent tubes. 
Leaf-discs from 84K were infected with Agrobacter- 
ium cultures containing Pptowusa-'GUS, Pptowox4a''GUS, 
P PtoWOXSa' 'GUSy Pptowoxii/i2a''GUSy 35S::PtoWUSa, 3SS:: 
PtoWOX4a, 3SS::PtoWOXSa, and 3SS:: PtoWOXll/12a 
constructs with ODeoo = 0.3-0.8. Infected leaf-discs were 
co-cultured with Agrobacteria in the AS induction 
medium [SIM; Murashige-Skoog (MS) basal medium with 
0.5 mg/1 6-benzyl aminopurine (6-BA) and 0.05 mg/1 
naphthaleneacetic acid (NAA)] in the dark for 3 days at 
22 ± 2°C. Leaf-discs were transferred on SIM with 3 mg/L 
hygromycin and 200 mg/L timentin under a 16 h/8 h 
(light/dark) regime. After one month, individual regener- 
ated shoots were removed and transferred onto root 
induction medium (RIM, 1/2 MS medium supplemented 
with 0.05 mg/L IBA, 0.02 mg/L NAA) containing 3 mg/L 
hygromycin and 200 mg/L timentin for AR induction. 
For AS induction, leaf discs were cultured on SIM for 
18 days (Additional file 1: Figure S6A-F), while stems 
were cultured on RIM for 9 days for AR regeneration 
(Additional file 1: Figure S6G-J). At least three independ- 
ently transformed lines were used for overexpression ana- 
lysis and promoter::GUS assay analysis. At least 25 clones 
of each overexpressed line were used for AR regeneration. 
All experiments were repeated at least three times 
with similar results. The regeneration of AS and AR 
in P. tomentosa was performed as above, but without 
the Agrobacteria co-cultured procedure. 

RNA isolation, RT-PCR, qRT-PCR, and RNA-seq analysis 

Total RNAs were extracted from roots, leaves, stem, 
CAM zone, and immature xylem of P. tomentosa Carr, 
and from materials at different stages and vegetative 
tissues of the hybrid poplar P. alba X P. glandulosa for 
RNA-sequencing using an RNeasy Plant Mini Kit and 
RNase-free DNase I set (Qiagen, Hilden, Germany). 
First-strand cDNA synthesis was carried out with 
approximately 1.5 (ig RNA using the Superscript III 
first-strand synthesis system (Life technologies, Carlsbad, 
California, U.S.) according to the manufacturer s instruc- 
tions. Specific RT-PCR primers were designed to have 
melting temperatures of 58-60°C and amplicon lengths 
of 150-260 bp using Primer3 software (http://frodo.wi. 
mit.edu/primer3/input.htm). The amplified fragments 
were separated on agarose gel electrophoresis. Real-time 
qRT-PCR was performed in quadruplicate using the SYBR 
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Premix Ex Taq™ II Kit (TaKaRa Dalian, Dalian, China) on a 
Roche lightCycler 480 (Roche Applied Science, Penzberg, 
Upper bavaria, Germany) according to the manufacturers 
instructions. Quantification was performed using Lightcy- 
cler 480 software (Version 1.5.1.62, Roche). Expression 
was normalized relative to the control (Actin) using 
Roche LightCycler advanced relative quantification analysis 
(E-method, http://www.roche-applied-science.com/shop/prod- 
ucts/gene-quantification-on-the-lightcycler-480-system) and 
fold changes (relatively to Bl or Rl stage) are shown in 
Figure 6. All experiments were repeated at least three times 
with similar results. 

PromoteniGUS assay 

Histochemical GUS staining was performed as follows: 
2-week old seedlings and 4-week old stem sections were 
first fixed in 90% cold acetone. Each sample was washed 
three times on ice using GUS staining buffer containing 
50 mM sodium phosphate (pH7.0), 2 mM potassium 
ferrocyanide, 2 mM potassium ferricyanide, 10 mM 
EDTA, and 0.2% (v/v) Triton X-100. Fixed samples were 
transferred into the staining solution [GUS staining buffer 
with 20% (v/v) methanol, and 1 mM X-Gluc] and slowly 
vacuumed. After 12 h incubation at 37°C with gentle 
agitation, samples were rinsed in 70% ethanol for imaging. 
At least five clones for every promoterr.GUS line were used 
for GUS staining. All experiments were repeated at least 
three times with similar results. 

Protein localization analysis 

Tobacco leaf samples were analyzed 48 to 72 h after 
infiltration. Fluorescence of plant epidermal cells was 
observed using a LSM 510 confocal laser scanning 
microscope (Carl Zeiss AG, Oberkochen, Germany). 
Detection of fluorescence was performed as follows: 
fluorescence imaging of YFP excitation at 514 nm; 
scanning at 520-555 nm. 

Availability of supporting data 

The phylogenetic trees generated in this study are available 
in Treebase under study number 15612 (http://treebase. 
org/treebase-web/search/study/summary.html?id= 15612) . 

Additional files 



genes in different tissues. A. Heat map of hierarchical clustering of PagWOX 
genes in vegetative tissues (YL, young leaves; ML, mature leaves; PS, primary 
stem; SS, secondary stem; R, roots). The data was obtained from our 
unpublished RNA-seq data. The expression level of genes was determined 
based on the value of RPKM (reads per kilobase of exon region in a gene 
per million mapped reads). The clustering was made on expression pattern. 
Details of the RPKM are shown in Table S3. Color scale represents log2 
expression values. Figure S4. The stages of AS and AR regeneration. A-F 
Leaf explants in indicated stages of AS regeneration. B1-B6 represent the 
status of explants on 0, 6, 9, 12, 15, 18 days after AS induction. G-J Stem 
segments in indicated stages of AR regeneration. R1-R4 represent the stem 
status on 0, 3, 6, 9 days after AR induction. Bar represents 1 cm in A-J. 

Additional file 2: Table SI. Primer sequences for amplification of 
PtoWOXs and promoters and qPCR analysis. Table S2. Similarity of WOX 
sequences from P. trichocarpa and P. tomentosa. Table S3. RPKM value of 
PagWOX genes in vegetative tissues obtained from RNA-seq data. Table S4. 
Accession number of proteins analyzed in this study. 
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